Coulomb dissociation of the deformed halo nucleus 31 Ne 
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The recently observed large cross sections for the Coulomb dissociation of 31 Ne nucleus indicate 
that this nucleus takes a halo structure in the ground state. We analyse these experimental data 
using the particle-rotor model, that takes into account the rotational excitation of the core nucleus 
30 Ne. We show that the experimental data can be well reproduced when the quadrupole deformation 
parameter of 30 Ne is around @2 = 0.2 ~ 0.3, at which the ground state takes the spin and parity 
of I n =3/2~. This state corresponds to the Nilsson level [330 1/2] in the adiabatic limit of the 
particle-rotor model. On the other hand, the state corresponding to the Nilsson level [321 3/2] with 
/?2 ~ 0.55 can be excluded when the finite excitation energy of the core is taken into account, even 
though this configuration is a candidate for the ground state of 31 Ne in the Nilsson model analysis. 
We discuss briefly also a possibility of the I 7r — l/2 + configuration with 02 ~ 1 and /?2 ~ —0.4. 

PACS numbers: 21.10.Gv,21.60.-n, 25.60.Gc, 27.30. +t 



A halo structure with an extended density distribu- 
tion is one of the characteristic features of weakly bound 
neutron-rich nuclei. It was first recognized in n Li, show- 
ing an abnormally large interaction cross section 
Other examples include 11 Be 0, Q and 19 C Q, both 
of which are typical one-neutron halo nuclei. The root- 
mean-square radius diverges for s and p waves as the 
single-particle energy approaches to zero [5[ , and the halo 
structure has been ascribed to an occupation of an I = 
or I = 1 orbit by the valence neutron [6[. 

The halo structure induces a large concentration of the 
electric dipole (El) strength in the low-excitation energy 
region, that is, a soft El excitation [1,0,9- Recently, the 
Coulomb breakup cross sections for 31 Ne were measured 
by Nakamura et a/.0], indicating a soft El excitation 
in the 31 Ne nucleus. Notice that a naive spherical shell 
model leads to the I/7/2 configuration for the valence 
neutron of 31 Ne. In order to generate the halo structure 
within the mean-field picture, the valence neutron there- 
fore needs to move in a deformed mean-field potential 
in which the s or p wave component in a weakly-bound 
single-particle wave function makes a dominant contribu- 
tion pL0TTl2l ] . Based on this idea, Hamamoto has carried 
out Nilsson model calculations with a deformed Woods- 
Saxon potential, and argued that the large Coulomb dis- 
sociation cross sections of 31 Ne can be understood if the 
valence neutron occupies [330 1/2], [321 3/2], or [200 1/2] 
Nilsson levels QJ] leading to F = 3/2" ([330 1/2] and 
[321 3/2] configurations) or 1/2+ ([200 1/2] configura- 
tion) for the spin and parity of the ground state of 31 Ne. 
We mention that the ground state with I" = 3/2~ is 
found also with shell model calculations [t| [13 as well as 
with a microscopic cluster model calculation }15|. 

In order to describe odd-mass deformed nuclei, espe- 
cially nuclei with small deformation as well as transi- 
tional nuclei, the particle-rotor model (PRM) has been 
applied with great success [l!14l8l |. It was also applied 
successfully to a halo nucleus 11 Be [H, [2(3] and neutron- 
rich Carbon isotopes 15 > 17 ' 19 C[2l|. In 30 Ne, the candi- 
dates for the first excited 2 + and 4 + states have been 



identified experimentally at excitation energies of 0.801 
MeV and 2.24 MeV, respectively [H, El- The energy 
ratio £4+ / E 2 +=2.80 suggests this nucleus to be a tran- 
sitional one in comparison with the ratio 3.33 for well- 
deformed nuclei. Although the experimental uncertainty 
is still large, an extremely small neutron separation en- 
ergy S„=0.29±1.64 MeVp4| suggests a halo structure of 
31 Ne. Because of these reasons, in this paper we assume 
a deformed 30 Ne core and adopt the particle-rotor model 
to study the ground state and the excited state properties 
of 31 Ne. 

The Nilsson model adopted in Ref. [lH to analyze 
the ground state configuration of 31 Ne corresponds to a 
strong coupling limit of the particle-rotor model. That 
is, the non-adiabatic effect due to the rotational states of 
the deformed core was neglected in the analysis. In gen- 
eral, in deformed nuclei with mass A< 40, the rotational 
energy is not small, and the non-adiabatic effect may play 
a role in discussing physical observables such as electric 
dipole (El) transitions. While the non-adiabatic effect 
may be estimated perturbatively in the Nilsson model, 
especially to determine the spin and parity of the nu- 
cleus in the laboratory frame [ill [l7j , it is not obvious 
whether the perturbation works for the particular case 
of 31 Ne nucleus. In this respect, it is highly desirable to 
disentangle how much the non-adiabatic effect appears 
in the assignment of the spin-parity of the ground state 
and also in the El transitions of 31 Ne. 

The aim of this paper is thus to study the 31 Ne nucleus 
by solving the particle-rotor model Hamiltonian without 
resorting to the adiabatic approximation and discuss the 
El excitation and the Coulomb breakup cross sections. 
In order to describe the structure of the 31 Ne nucleus, 
we assume that the core nucleus 30 Ne is statically de- 
formed with a quadrupole deformation parameter 02 and 
the motion of the valence neutron is coupled to the rota- 
tional motion of the core nucleus. Assuming the axially 
symmetric shape of the core nucleus, we consider the fol- 
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lowing Hamiltonian for the ciu Ne+n system: 

h 2 



H 



2/i 



V" + V(r, f c ) + H IC 



(1) 



where fi — m-^A c /{A c + 1) is the reduced mass of the 
valence neutron with A c —30 and win being the mass 
number of the core nucleus and the nucleon mass, re- 
spectively. -ff r ot is the rotational Hamiltonian for the 
core nucleus given by H mt = I 2 h 2 /2J, where I c is the 
spin of the core nucleus and J is the moment of inertia. 
V is the single-particle potential for the valence neutron 
interacting with the deformed core, r and r c are the 
coordinate of the valence neutron and the direction of 
the symmetry axis of the core nucleus in the laboratory 
frame, respectively. For the single-particle potential V, 
we use a deformed Woods-Saxon potential, 



where / is the spin of 3 Ne and M is its z-component. 
The expansion coefficients as weu as ^ ne corre ~ 

sponding eigenenergy for the 31 Ne nucleus are obtained 
by numerically diagonalizing the Hamiltonian H. 

In order to identify the ground state configuration, we 
first solve the Hamiltonian by setting H mt — in Eq. 
([T]). In this case, the K quantum number, that is, the 
projection of the total angular momentum onto the z-axis 
in the body-fixed frame, is conserved, and several states 
with different /, having the same value of K, are degen- 
erate in energy when the maximum value of I c included 
in the calculation is sufficiently large. As has been shown 
in Ref. [2(|, the wave function in this limit is related to 
the wave function in the Nilsson model, 4'jiK, as, 



V(r,r c ) = -V ws [ l-F ls r 2 (l-s)-4- ) f{r) 



r dr 



df(r) 

+V ws Rot3 2 J Y 1 Y 20 (r cn ) 7 
dr 

V (r) + V 2 (r)Y 20 (r CD ), 



(7) 



where 
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2 ' ■ + ±.V2(jKI c 0\IK), (8) 
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where f cn is the angle between r and f c , and f(r) is given and 4'jiK satisfies [27 
by 



/(r) = l/[l + exp((r-i2o)/o)] 

,1/3 



(4) 



with i?o = tqA c ■ Here, for simplicity, we have ex- 
panded the deformed Woods-Saxon potential up to the 
linear order of /3 2 . We have checked that this approx- 
imation works for the range of deformation parameter 
considered in this paper, by comparing our results for 
the Nilsson model with the numerical results obtained 
with the code WSBETAp^. 

In principle, the eigenfunctions of the Hamiltonian ([1]) 
can be obtained by solving the coupled-channels equa- 
tions in the coordinate space both for the bound and the 
scattering states [13, |2lj, [26|. A continuum B(E1) dis- 
tribution can then be constructed with those scattering 
wave functions. However, since the experimental data 
are so far available only in the form of inclusive cross 
sections, we will not discuss the detailed structure of the 
strength distribution here but leave it for future publica- 
tions. It is sufficient for our purpose to expand the wave 
function on some basis and compute a discrete strength 
distribution. We do this with the eigenfunctions of the 
spherical part of the potential V, 
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V(r, 
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•jiK{r)yjiK{r) 



(9) 

One can regard Eq. ([7]) as a transformation of the Nilsson 
wave function from the body-fixed frame to the labora- 
tory frame, where the angular momentum is conserved. 
Notice that 
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2/x 



V 2 + VoM -e nji 



RnjiirWjimir) = 0, (5) 



where R n ji(r) is the radial wave function and yji m (r) 
is the spin-angular function. The continuum spectrum 
can be discretized within a large box. Together with 
the rotational wave function 0/ c ^/ c (f c ), the total wave 
function for the n+ 30 Ne system is expanded as, 

*/A/(r,r c ) = ^^«S /c i?^M[^(r)0/ c (f c )]( /M ), 

njl I c 

(6) 



is the probability of each (J, I) component in the Nilsson 
wave function and is independent of / when H Iot = 0. 

In this limit, we therefore obtain a collection of Nilsson 
levels. As usual, we put two neutrons to each Nilsson 
orbit from the bottom of the potential well, and seek 
the Nilsson orbit which is occupied by the last unpaired 
neutron. We then gradually increase the value of the 
2 + energy of the core nucleus up to the physical value, 
E 2 + =0.801 MeV, and monitor how the Nilsson orbit for 
the valence neutron evolves. For a finite value of E 2 + , the 
K quantum number is not conserved any more due to the 
Coriolis coupling, and the degeneracy with respect to I is 
resolved. We select the lowest energy state among several 
/ at E 2 + =0.801 MeV as the ground state of 31 Ne. In this 
way, we take into account the Pauli principle between the 
valence neutron and the neutrons in the core nucleus. 

Table I shows the probability of each component, 

PjUc = J2n\ a njii c \ 2 ' in tne wave function for ,02=0.1, 
0.2 and 0.55 with two different values of the ground state 
energy of 31 Ne. The ground state energy obtained by 
diagonalizing the Hamiltonian H is measured from the 
threshold of n+ 30 Ne, and thus it is equivalent (except 
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TABLE I: The probability of the component [I c ® (jl)] in deformed wave functions of 31 Ne, where I c is the spin of the rotational 
state of the core nucleus 30 Ne and (jl) is the angular momentum between the valence neutron and the core nucleus. These 
are obtained with the particle-rotor model (PRM), where the results with vanishing rotational energy, E 2 +=0, correspond 
to the Nilsson model calculations. The states with ,32=0.1 and 0.2 correspond to the Nilsson orbit [330 1/2] while the state 
with p 2 =0.55 corresponds to the Nilsson orbit [321 3/2]. The K 7 ' quantum number is therefore K 7r — l/2~ and K n =3/2~ , 
respectively. The spin of 31 Ne in the laboratory frame is I n = 3/2" for f3 2 =0.2 and 0.55 while it is 7/2" for /3 2 —0.l. Two 
different values of the one neutron separation energy S n of the 31 Ne nucleus are compared. For the case of E 2 + = 0, we also 
list the total p 3 / 2 and f 7 / 2 probabilities defined as Pji = J2i Pju c - 





[0+® p 3 / 2 ] 


[2 + ® P3/2] 


(total p 3/2 ) 


[2+® P1/2] 


[0+<S> f 7/2 ] 


[2+® f 7/2 ] 


[4+® f 7/2 ] 


(total f 7/2 ) 


S n 


=0.3 MeV 






















=0.1, E 2+ 


= (Nilsson) 


0% 


19.5% 


30.3% 


0% 


17.2 % 


20.5% 


18.1 % 


68.8% 




E 2 + 


= 0.801 MeV (PRM) 


0% 


6.5% 




% 


75.2% 


17.5 % 


0.5 % 






=0.2, E 2+ 


= (Nilsson) 


22.2% 


22.2% 


44.4% 


2.9% 


0% 


33.5 % 


18.6 % 


52.1% 




E 2 + 


= 0.801 MeV (PRM) 


44.9% 


8.4% 




2.0% 


0% 


42.7 % 


1.5 % 




fa-- 


=0.55, E 2 + 


= (Nilsson) 


11.4% 


11.5% 


22.9% 


0% 


0% 


25.5 % 


45.9 % 


71.4 % 




E 2+ 


= 0.801 MeV (PRM) 


1.9% 


29.7% 




4.4% 


0% 


23.0 % 


35.5 % 




S n 


=0.1 MeV 




















fa 


=0.1, E 2+ 


= (Nilsson) 


0% 


26.7% 


41.5% 


0% 


14.3% 


17.0% 


15.0% 


57.1% 




E 2 + 


= 0.801 MeV (PRM) 


0% 


7.2% 




0% 


74.8% 


17.1% 


0.5% 




/?2 = 


=0.2, E 2+ 


= (Nilsson) 


26.5% 


26.5% 


53.1% 


3.8% 


0% 


27.4 % 


15.2 % 


42.6% 




E 2 + 


= 0.801 MeV (PRM) 


53.4% 


8.3% 




2.1% 


0% 


34.6 % 


1.2 % 




ft = 


=0.55, E 2+ 


= (Nilsson) 


12.9% 


12.9% 


25.7% 


0% 


0% 


24.6 % 


44.2 % 


68.8 % 




E 2 + 


= 0.801 MeV (PRM) 


1.8% 


30.9% 




4.9% 


0% 


23.6 % 


33.5 % 





for the negative sign) to the one-neutron separation en- 
ergy, S n (it should not be confused with the energy e n ji in 
Eq. |S1)). Following Ref. 0, we use r =1.27 fm, a=0.67 
fm, and i<i s =0.44 for the parameters of the Woods-Saxon 
potential, V\ys is adjusted to reproduce a given value 
of S n . We use i?box=60 fm for the size of the box to dis- 
cretize the continuum spectrum, and include the single- 
particle orbits © up to e = 90 MeV and Z max =10. For 
the spin of the rotational states of the core nucleus, we 
include up to I c = 8. In the case of vanishing rota- 
tional energy, the states with /3 2 =0.1 and 0.2 correspond 
to the Nilsson orbit [330 1/2] {K* = 1/2") while the 
state with ,$2=0.55 corresponds to the Nilsson orbit [321 
3/2] (K n = 3/2"). See e.g., Fig. 2 in Ref. [H for 
the Nilsson diagram. The spin of 31 Ne in the labora- 
tory frame is T* = 3/2" for /? 2 =0.2 and 0.55, while it is 
r '= 7/2" for /3 2 =0.1. 

In the adiabatic limit (that is, the Nilsson model), 
for a given j and I, the contribution of each compo- 
nent of [I c <g) (jl)] with different value of I c to the El 
excitation is similar to each other since the radial de- 
pendence of the wave function is identical (see Eq. (0). 
Therefore, the relevant quantity in this limit is the to- 
tal probability for each j and I given by Eq. (Ti"Uj) . On 
the other hand, in the case of finite rotational energy 
(i.e., the non-adiabatic coupling), the radial dependence 
of the wave function largely depends on I c . The wave 
function is spatially most extended for the component 



which couples to I c = 0, as the absolute value of the 
diagonal energy E — E Iot is the smallest (that is, most 
weakly bound). Therefore, the relevant quantity in this 
case to the El excitation is the probability of the I c = 
component. In the case of /3 2 = 0.2, the probability of 
[0 + (g>p 3 / 2 ] increases considerably as the rotational energy 
increases from zero to 0.801 MeV, and the probability at 
E 2 + = 0.801 MeV is similar to the total p 3 / 2 probability 
in the adiabatic limit. As the separation energy decreases 
from S"„=0.3 MeV to SW^O.l MeV, the total n-wave prob- 
ability increases in the adiabatic limit |10l - [l3l ]. and cor- 
respondingly the probability of [O+Cgipa^] also increases 
in the non-adiabatic case. In contrast, in the case of 
(3 C = 0.55, the probability of [0 + £g>p 3 / 2 ] decreases consid- 
erably as the rotational energy increases. The probabil- 
ity remains small even if the separation energy decreases 
from S" n =0.3 MeV to S„=0.1 MeV. We have confirmed 
that this is the case even when the separation energy is as 
small as S n = 0.01 MeV. We thus expect that the disso- 
ciation cross section decreases significantly for this con- 
figuration due to the non-adiabatic effect. In the case of 
/3 2 = 0.1, the [0 + <8>f 7 / 2 ] component increases significantly 
for the non-adiabatic coupling. Since the /-wave does not 
form a halo structure, the El excitation probability will 
decrease considerably if the rotational energy is taken 
into account. As the separation energy decreases, the p- 
wave component will increase in the Nilsson wave func- 
tion. However, the non-adiabatic effect always quenches 
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FIG. 1: (Color online) The Coulomb dissociation cross sec- 
tions of the 31 Ne nucleus as a function of the one- neutron 
separation energy S n . The upper panel shows the results of 
the Nilsson model calculations, while the lower panel shows 
the results with a finite rotational energy of the core nucleus. 
The dotted, the solid, and the dashed lines denote the results 
with fa — 0.1, 0.2, and 0.55, respectively, where the spin and 
the parity of the ground state is 7/2~, 3/2~, and 3/2~, re- 
spectively. The shaded region indicates the experimental data 



the p-wave component in this weak coupling regime, and 
one could not expect large Coulomb dissociation cross 
sections with this configuration. 

Let us now numerically compute the Coulomb disso- 
ciation cross sections and confirm the above behaviors 
discussed in the previous paragraph. To this end, we 
first compute the El strength, 



B(El;i->f) 



1 



2I t + 1 



D 



(11) 



where the initial and final wave functions, ^ and 'F/, 
have the same form as in Eq. (j6|) , and the dipole operator 
is given by = -[Z c e/(A C + 1)] • rY^(r). The 
Coulomb dissociation cross sections can be obtained by 
multiplying the El virtual photon number, Nei, to the 
B(E1) strength [281-130}. Summing all the final states, the 
inclusive Coulomb breakup cross section reads, 



16tt 3 
9Hc 



N m (E f -Ei)-B(El;i->f). (12) 



Since we compare our results with the experimental data 
for the dissociation cross sections, we restrict the summa- 
tion in Eq. (|T2]) only to those states above the threshold. 

Figure 1 shows the Coulomb dissociation cross sections 
for the three configurations listed in Table I. Since the 
empirical separation energy S n has a large uncertainty, 
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FIG. 2: (Color online) The Coulomb dissociation cross sec- 
tions of the 31 Ne nucleus as a function of the one-neutron 
separation energy S„ at /32=0.95. The corresponding Nils- 
son orbit is [200 1/2] (K n = 1/2+), and the spin-parity of the 
ground state is I n — l/2 + . The dashed and the solid lines de- 
note the results in the adiabatic limit and in the non-adiabatic 
case, respectively. 



we plot the calculated cross sections as a function of S n . 
The upper panel shows the results obtained in the adia- 
batic limit, while the lower panel shows the results with 
the finite rotational energy. The dotted, solid, and the 
dashed lines are obtained with fa— 0.1, 0.2, and 0.55, re- 
spectively. The shaded region indicates the experimental 
data @. One can clearly see in the lower panel of Fig. 1 
that the experimental data can be well reproduced with 
fa = 0.2 if the separation energy S n is around 0.17 MeV. 
The results remain almost the same even if we vary the 
deformation parameter in the range of 0.17 < fa < 0.33. 
On the other hand, the cross sections obtained with fa= 
0.1 and 0.55 appear too small to account for the experi- 
mental data. These configurations may yield a reasonable 
reproduction of the experimental data in the adiabatic 
limit because of the p-wave dominance in the wave func- 
tion, but it is simply an artifact of the adiabatic approxi- 
mation. The non-adiabatic effect eliminates the possibil- 
ity of these configurations to be the ground state of 31 Ne. 
All of these behaviors agree with the expectations. 

Besides the configuration with fa ~0.2, there may be 
other possibilities for the ground state configuration of 
31 Ne. It was argued in Ref. [H that the r = 1/2+ 
configuration originated from the [200 1/2] Nilsson or- 
bit {K* = l/2 + ) cannot be excluded if the deformation 
parameter fa is large. Figure 2 shows the Coulomb dis- 
sociation cross sections calculated with fa = 0.95. The 
dashed and the solid lines correspond to the results in 
the adiabatic limit and those with the finite rotational 
energy, respectively. For this configuration, the cross 
sections show some jumps when a resonance state be- 
comes a bound state as the separation energy increases 
so that its contribution is excluded in the sum in Eq. 
(fT2)) . In any case, the experimental cross sections can 
be reproduced when the separation energy is around 0.8 
MeV, and thus this configuration cannot be excluded, al- 
though it is required that 31 Ne has a bound excited state 



5 




FIG. 3: (Color online) Same as Fig. 2, but for /3 2 = -0.35. 
The corresponding Nilsson orbit is [200 1/2] (K n = 1/2+ ), 
and the spin-parity of the ground state is T* = l/2 + . 



and/or low-lying resonance states. For this configuration, 
the [0 + <g)Si/2] configuration is the main component of the 
wave function both in the adiabatic limit and in the non- 
adiabatic case (for instance, the probability is 80.5 % in 
the adiabatic limit and 79.2 % in the non- adiabatic case), 
and the non-adiabatic effect on the Coulomb dissociation 
is found to be small. 

Although it may be unlikely that the 30 Ne nucleus is 
oblately deformed [23| , the Coulomb dissociation data by 
themselves do not exclude a possibility of the oblate con- 
figuration corresponding to the [200 1/2] {K* = 1/2+) 
Nilsson diagram at fa ~ —0.35. In order to demonstrate 
this, Fig. 3 shows the calculated Coulomb dissociation 
cross sections obtained with fa = —0.35. The spin-parity 
of the ground state is T" = l/2 + . As one can see, the 
experimental Coulomb dissociation cross sections can be 
well reproduced when S n ~ 0.2 MeV. The main compo- 
nent of the wave function is the [2 + £g)d3/2] configuration 



with an appreciable mixture of the [0 + Cg)Si/2] configura- 
tion. The s-wave probability does not change much by 
the non-adiabatic effect, and the experimental data can 
be reproduced both in the adiabatic limit and in the non- 
adiabatic case. 

In summary, we have discussed the El excitation of the 
31 Ne nucleus using the particle-rotor model. The finite 
rotational excitation energy of the core nucleus has been 
taken into account. We have shown that the experimen- 
tal cross sections can be well reproduced with the defor- 
mation parameter of 0.17 < fa < 0.33 and the separation 
energy of 0.13 < S n < 0.2 MeV, for which the ground 
state configuration has the spin-parity of I" = 3/2~. 
However, the I n = 1/2+ configurations with large de- 
formation around fa ~ 0.95 or with oblate deformation 
around fa ~ — 0.35 cannot be excluded only from the 
Coulomb dissociation data. On the other hand, we have 
shown that the J* = 7/2 ~~ configuration at fa ~ 0.1 and 
the -P = 3/2 ~ configuration at fa ~ 0.55 can be excluded 
due to the non-adiabatic effect, even though those con- 
figurations may reproduce the experimental data in the 
adiabatic limit. 

For the 31 Ne nucleus, the momentum distribution for 
the nuclear breakup process [3l[ as well as the interac- 
tion cross section have been recently measured at the 
RIBF facility at RIKEN [Hill]. It will be interesting to 
investigate whether our model can reproduce these ex- 
perimental data simultaneously. We will report on it in 
a separate publication. 
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